Plot[97.6215 Sqrt[ r^2 (BesselI[0, 0.185185 r] BesselK[0, 0.185185 r] - BesselI[1, 0.185185 r] BesselK[1, 0.185185 r])]], {r, 0, 40}]